Persistence in g-state Potts model: A Mean-Field approach 
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We study the Persistence properties of the T = coarsening dynamics of one dimensional q- 
state Potts model using a modified mean- field approximation (MMFA). In this approximation, the 
spatial correlations between the interfaces separating spins with different Potts states is ignored, 
but the correct time dependence of the mean density P(t) of persistent spins is imposed. For this 
model, it is known that P(t) follows a power-law decay with time, P(t) ~ t~°^ where 9{q) is the 
g-dependent persistence exponent. We study the spatial structure of the persistent region within 
the MMFA. We show that the persistent site pair correlation function P2{r,t) has the scaling form 
P2(r,t) — P(t) 2 f(r /tz ) for all values of the persistence exponent 9(g). The scaling function has 



the limiting behaviour f(x) 



(x <C 1) and f(x) — > 1 (x ^> 1). We then show within the 



Independent Interval Approximation (IIA) that the distribution n(k, t) of separation k between 
two consecutive persistent spins at time t has the asymptotic scaling form n(k,t) — t~ 2< ^ 'git, i) 
where the dynamical exponent has the form 0=max(i,6>). The behaviour of the scaling function 
for large and small values of the arguments is found analytically. We find that for small separations 
k <C v ,n(k,t) ~ P(t)k~ T where r=max(2(l — 9), 29), while for large separations k ^> r, g{t,x) 
decays exponentially with x. The unusual dynamical scaling form and the behaviour of the scaling 
function is supported by numerical simulations. 



I. INTRODUCTION 



In recent times, the notions of Persistence and the asso- 
ciated Persistence exponent has become one of the active 
topics of research in Non-equilibrium Physics [1] . In gen- 
eral, the persistence probability P(t) is the probability 
that a stochastic variable <f>{t) remains above or below a 
certain arbitrary value (say, its initial value) for a time 
interval [0 : t] . The idea of Persistence is particularly rel- 
evant in coarsening systems where P(t) has the natural 
interpretation of being the fraction of space in the system 
which remains in the same phase, starting from a random 
initial configuration. In these systems (as in some other 
systems also) P(t) typically decays as a power-law at 
large t, with an exponent that is non-trivial to compute 
and does not appear to be simply related to other known 
exponents that characterize the process. 

The coarsening dynamics of the one dimensional q- 
state Potts model at zero temperature is one of the few 
cases where the Persistence exponent 9(q) is known ex- 
actly. The solution was provided by Derrida et. al through 
a mapping of the process to a coagulation model in steady 
state [2]. It was shown that at late times t, well beyond 
the time scale of equilibration, the fraction of persistent 
spins left in a finite system of linear dimension L scale 
as P q (t > L 2 ,L) ~ L~ 2e( - q \ where 6(q) is given by the 
non-trivial expression 



%) 
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(1) 



For times t < L 2 , it follows that P q (t) ~ t^ 9 ^. An 
interesting question in this context is about the spatial 
distribution of spins which are persistent up to any given 
time t. Clearly, in a many-body process like the time 
evolution of Potts model, the probability that a given 
spin is persistent is closely linked to the state of other 
spins. This inter-dependence of spins is crucial in that it 
makes the time evolution at any single site strongly non- 
Markovian, which makes the computation of the Persis- 
tence exponent highly non-trivial. A study of the spatial 
aspects of the Persistence problem is therefore impor- 
tant from the point of view developing an intuitive un- 
derstanding of the phenomenon, and also illuminates the 
interplay between Persistence decay and the underlying 
domain coarsening process. 

The spatial distribution of persistent sites and its time 
evolution have been studied previously through numeri- 
cal simulations in one dimensional diffusion equation [3] , 
q-state Potts model [4-7] , two dimensional Ising model 
[8] and one dimensional Ising model with parallel dynam- 
ics [9] . An analytic study using a rate equation approach 
under the Independent Interval Approximation(IIA) has 
been carried out for one dimensional A + A — > model 
which is closely related to ID Ising model [6]. It is 
now generally understood from physical arguments and 
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simulations that for a coarsening process in d dimen- 
sions where characteristic length scale has the power-law 
growth L{t) ~ t 1 / 2 , the set of persistent sites forms a 
fractal structure with fractal dimension df ~ d — z6 over 
length scales r <C t 1 ^ [7,10]. The distribution is ho- 
mogeneous beyond this length scale. Furthermore, since 
df > 0, the distribution is expected to be homogeneous 
over all length scales if 9 > d/z. This has important 
consequences for systems like the Potts model where 9 
changes with the Potts state q. In particular, for Potts 
model in d = 1, Bray and O'Donoghue [7] argued that a 
transition from fractal to homogeneous distribution oc- 
curs as 8 crosses ^ ■ This transition is also marked by an 
abrupt change in the dynamical exponent characterizing 
the separation between persistent sites. The character- 
istic length scale was conjectured to have the unusual 
dynamical scaling form C(t) <~ with <p = max(^,9). 
This conjecture based on physical arguments was sup- 
ported by numerical simulations [7] . 

In this paper, we use a Mean-Field approach to address 
the problem of spatial distribution of persistent sites in 
g-state Potts model. The essential idea behind this ap- 
proach is as follows. It is well-known that the T = 
coarsening dynamics of the g-state Potts model can be 
mapped to a reaction-diffusion process. In this process, 
the interfaces between different species of Potts spins 
are represented by diffusing particles A, which annihilate 
or coagulate upon meeting with g-dependent probabili- 
ties. In the mean field approach, these diffusing particles 
are treated as homogeneously distributed, with (time- 
dependent) density equal to the average density in the 
original reaction-diffusion problem. This approach has 
been discussed in some earlier works as a heuristic argu- 
ment [11] and as a toy model for persistence [12]. We 
argue that this approach yields a lower bound for the 
Persistence exponent in the Potts model. We then con- 
struct an artificial model which is devoid of spatial cor- 
relations among diffusing particles, but with persistence 
exponent tuned to be exactly equal to the Potts model 
value. We refer to this model as the Modified Mean-Field 
Approximation (MMFA) and use this approximation to 
study the spatial distribution of persistent sites in g-state 
Potts model. 

We outline our main results at this point. Within the 
MMFA, we show analytically that the correlation length 
for the spatial distribution of persistent spins scales as 
~ fa and the equal time pair correlation function 
Pi (r, t) (defined as the probability that the spin at ori- 
gin and the point r are persistent at time t) has the 
scaling form P 2 (r, t) = P (t) 2 / (r / \fi) for any value of 
9. This shows that the persistent spins have a fractal 
distribution with df = 1 — 29 over length scales r < t5 
when 9 < \, but df = when 6 > |. We find that the 
characteristic length scale of the spatial distribution of 
persistent spins has the unusual scaling form £(t) ~ 
where </>=max(|, 0), in agreement with the conjecture in 



[7] . The empty interval distribution itself has the scaling 
form n(k, t) = C(t)- 2 g[t, k/C(t)], where g(t, x) ~ t-^x~ T 
for x <C 1 and decays exponentially with x when x » 1. 
The exponents ip and r depends on 9 through the rela- 
tions i> = 9{29 - 1)H(9 - i) and r=max[2(l - 6), 26], 
where H(x) is the Heaviside step function. We support 
these results with numerical simulations. 

The paper is arranged as follows. In the next section, 
we outline the mean field approach. In Section III, we in- 
troduce the MMFA and compute the pair correlation and 
Empty Interval Distribution of persistent sites to charac- 
terize their spatial distribution. These predictions arc 
compared with the results of numerical simulations in 
the (/-state Potts model in Sec. IV. We summarize our 
results and present our conclusions in Sec. V. 



II. THE MEAN-FIELD APPROXIMATION 

In the zero temperature coarsening dynamics of g-state 
Potts model in d = 1, the interfacial points between dif- 
ferent species of Potts spins perform independent ran- 
dom walks on the lattice and annihilate each other with 
probability )i or coagulate with probability ■ 

In the process, persistent sites are 'wiped out', and the 
surviving random walkers build up spatial correlations 
among themselves. The distribution of intervals between 
the surviving random walkers at any (sufficiently late) 
time t is described by a (stationary) scaling function 
which is known exactly for all values of q [13]. The aver- 
age density n(t) at time t decays as 



n(t) 
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(2) 



asymptotically [13]. The essential idea behind the 
mean field (MF) approximation is to treat the random 
walkers as forming a homogeneous background of average 
density n(t), as far as the persistent sites are concerned. 
We define the Persistence probability P(t) as the proba- 
bility that the site at the origin is unvisited by any walker 
till time t. Then, the probability that the site the prob- 
ability that the site is visited by a walker between time t 
and t + dt is — d . Let R{x, t) be the probability that 
the site at origin is visited by a walker for the first time at 
time t, whose initial position was x at t = 0. Within the 
MF approximation, any walker would survive with prob- 
ability n(t), and the probability that it will make a first 

visit to origin at time t is given by q(x, t) = \J^j§/2 e ~^*' 
[14]. It follows that R(x,t) = n(t)q(x,t). We now in- 
tegrate R(x,t) over all initial positions x and multiply 
by the probability that the origin is persistent at time t, 
which is simply P{t). So we find 



dP(t) 
dt 



= -P^n^K^t) 



d=l 



(3) 
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where K\(t) = q(x, i)dx is the Smoluchowski con- 
stant [15] in d = 1. After substituting for q(x, t) and n(t), 
we find = -%-P(t), so that P(t) ~ t-°'^\ and 

#*(<?) = -^f^f- is the Persistence exponent within the 
MF model [12]. It is interesting to compare the mean 
field prediction for 9 with the exact value of the expo- 
nent. For q = 2, 9* (2) ~ 0.225, while the exact value 
from Eq.l is |. For the q — oo case, the MF model pre- 
dicts 9*(oo) ~ 0.45, which is to be compared with the 
exact value 9(oo) = 1. Upon extending the comparison 
to the entire range of values of q, it is clear that the mean 
field treatment consistently under-estimates the value of 
9. 

We now argue that 9*(q) is, in fact, a lower bound 
for 9{q). In the mean field approach discussed so far, it 
is assumed that the random walkers disappear from the 
lattice at random at such a rate so that their average 
density falls as n(t). The actual reaction-diffusion pro- 
cess is quite different, because only walkers which come 
very close to another walker are likely to be removed. 
Clearly in regions of space where walkers come close to 
each other, they are likely to visit the same site again 
and again. This effect is much more within the mean 
field approach, where the walkers actually pass through 
each other, possibly several times before disappearing. 
Thus, it is plausible that for the same average density of 
walkers, a larger number of persistent sites will be vis- 
ited in the actual reaction-diffusion model, compared to 
its mean field analogue. Since this is true for all times, 
the average density of persistent sites in mean field theory 
will be higher than the same in the actual Potts model 
dynamics. This would naturally imply that 



0*(q) < 9{q) 



(4) 



Interestingly, we show now that the Mean Field ar- 
gument presented above yields the correct value for the 
Persistence exponent for A + A — > model in d = 2. This 
is not surprising, since for this model, the upper critical 
dimension is d c = 2, and the mean field treatment be- 
comes exact above this dimension. It can be shown that 
in d = 2, the probability that a random walker starting at 
an an arbitrary point crosses a circle of radius a around 
the origin for the first time at t is given by the expression 
[16] 



logiADtja 2 ) 



t > a 2 /D 



(5) 



which is the analogue of Smoluchowski constant in d = 
2. The asymptotic particle density decay for A + A — > 
model in d = 2 has the form n(t) ~ gj^ t^]' Upon 
extending the MF arguments presented previously, we 
find that 



After substituting for n(t) and K%(t), and taking the 
limit a -> 0, we find that P(t) ~ t~? so that 6* = 1/2 
in d = 2. This result is exact, as has been shown by a 
rigorous field-theoretic calculation [11]. 

We thus observe that while the mean field approach, in 
general, gives only a lower bound for the persistence ex- 
ponent, it correctly identifies the essential features that 
brings about this power-law decay, ie., the diffusive mo- 
tion of interfacial points and the decay in their over- 
all density. In the following sections, we use a slightly 
modified version of this treatment to study the spatial 
distribution of persistent sites in g-state Potts model. 



III. THE MODIFIED MEAN FIELD 
APPROXIMATION (MMFA) 

Our purpose is to use the mean field approach to study 
the spatial distribution of persistent spins in the g-state 
Potts model, and in particular, to understand the transi- 
tion from fractal to homogeneous distribution as 9 crosses 
i. However, it may be noted that in the mean field 
approximation to the dynamics of the Potts model, the 
largest value of 9 (attained at q = oo) is ~ 0.45 which 
is less than the transition value \. This problem is cir- 
cumvented by defining an artificial model where we define 
the diffusing particles as non-interacting random walkers, 
who can pass through each other. The model also allows 
for multiple occupancy of lattice sites. The dynamics 
consists of random walkers being picked at random and 
taken out of the lattice at a time-dependent rate, which 
is tuned to produce power-law decay P(t) <~ t~ e with 
any arbitrary 9'. From the arguments presented in the 
previous section, it is obvious that this would be the case 
if the average density were to decay as n(t) ~ -\/¥v? 
asymptotically. By construction, this model is devoid 
of spatial correlations among reacting particles (ie., it is 
still mean field) but 9' is now arbitrary. If we now choose 
9' = 9(q), this model is an approximation to the g-state 
Potts model, with the simplifying feature that the spa- 
tial correlation between interfacial points is now absent. 
We shall refer to this model as the Modified Mean Field 
Approximation (MMFA) for the original Potts model. 



dPjt) 
dt 



= -P(t)n(t)K 2 (t) 



d = 2 



(6) 



A. Pair correlation for persistent spins within the 
MMFA 



In this section, we compute the equal time pair corre- 
lation function for persistent spins in g-state Potts model 
under the MMFA. We define P2(r,t) as the probability 
that both the site at origin and the site at r > are 
persistent at time t. Our purpose is to compute P%{r, t) 
for various r. 

The generalization of Eq.3 to this case is 
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dP 2 (r,t) 
dt 



2P 2 (r,t) f 

J — ( 



R r (x,t)dx 



(7) 



where R r (x,t) is the probability that a particle with 
initial position x(— oo < x < r) will make a first visit 
to the origin at time t, without ever crossing r in the 
interval [0:tj. The factor 2 in front takes into account 
the probability that either of the sites could be reached 
by one of the diffusing particles. Unlike the first case, 
R r (x,t) is now different for x < and < x < r. For 
x < 0, the constraint of no crossing at r is irrelevant for 
the computation of R r (x, t) , since to reach r, the particle 
would have to cross the origin first. So R r (x, t) — R(x, t) 
simply for x < 0, and so 



J — ( 



R r (x, t)dx = — 
zt 



(8) 



For x > 0, this is no longer true, and R r {x, t) needs to 
be computed separately. The quantity that we need here 
is q r (x,t), the probability that a diffusing particle whose 
position at t = is x, will reach the origin for the first 
time at t, without ever crossing the point r in between. 
Then R r (x,t) — n(t)q r (x,t). To find q r (x,t), let us use 
the following standard method. Consider a random walk 
starting from 0<x<rati = with absorbing barriers 
at and r. If the probability distribution of the position 
z of the walker at time t is u x (z, t), then 



du x (z,t) 

Qr{x,t) = — | I=C 



(9) 



The expression for u x (z,t) is known exactly, and the 
asymptotic form at large t [14] is 



c (z,t) = 



E 



— e 



(10) 



with n — e X 2 . 
in Eq. 7, we find 



After substitution of Eq.8 and Eq.12 




FIG. 1. Plot of G(x) vs x 



dP 2 (r,t) 
dt 



-m{r,t) 



2t 



which admits a scaling solution of the form 



P 2 (r,t)=P(tff(-=) 



(14) 



(15) 



and the scaling function f(x) is given by the following 
expression. 



from which, we find 



q r (x,t) 



k— — oo 



X + 2kr (z-x-2fcr) 2 

e 2t 

t 



x — 2kr (z+x-2kT) 2 

1 e (11) 



We note that for r 3> 1 2 , the k = mode is the 
dominant term in the sum, and this gives q r (x,t) 



^p-xe 2t + smaller terms that vanish as — * °°- 
Clearly in this limit, we recover the r = 00 term, as we 
should. It then follows that 

J\ {X ,t) d x = l t G^) (12) 

where 



■2k) 2 



(13) 



2 dx 



(16) 



Let us now consider the limiting behaviour of the scal- 
ing function for x <C 1 and x >• 1. In the first case, it is 
clear from Fig.l that G(x) « 0, and so ||| = —9f(x), 
which implies that f(x) ~ x~ 2e as x — > 0. In the op- 
posite extreme G(x) — > 1 as x — > 00, and so |^ w 0, 
which means that /(x) approaches a constant value in 
this limit. It is further clear that, from the definition 
of the scaling function as given by the expression Eq.15, 
this constant is unity, since we expect P 2 (r,t) — > P(t) 2 
as r — > 00. For convenience of later calculations, we ap- 
proximate the scaling function as 



29-23 



f(x) - a M x 

m = 1 



: x < a 
: x > a 



(17) 



fe=i 



where a is a number, of order unity. 



4 



We see that under the MMFA, the pair correlation 
function has a scaling form which is same for all values 
of 9, with power-law decay P 2 (r,t) ~ P{t)r~ 26 for short 
distances r«(5. As is well-known, power law decay of 
pair correlation function points to the underlying scale in- 
variance of the spatial distribution of the persistent spins. 
This is characteristic of a fractal distribution under some 
circumstances. To see this, let us first define C(r, t) = 
P(t)~ 1 P 2 (r,t), which is the probability of finding a per- 
sistent spin at a distance r from another persistent spin. 
Now, the integral M(R,t) = C(r,t)dr is the total 
number of persistent spins within a radius R of a persis- 
tent spin. Clearly, from the scaling form described above, 
M(R,t) ~ R d ! for R < ti, where d f =max(l - 29,0). 
For R > ti, we find that M(R,t) ~ RP(t), which is 
simply a homogeneous distribution. Thus, if we look at 
length scales i? < ^, there is a fractal structure when 
9 < \. However, when 9 > |, this scale- invariant struc- 
ture is replaced by a few isolated sites, whose number 
does not grow with the length scale of observation. 

Clearly, the spatial distribution of persistent spins un- 
dergoes a transition as 9 crosses \. Indeed, if we consider 
time scales beyond equilibration time t » L 2 , for 9 < \ 
the total number of persistent spins left in the system 
scales as L 1 ^ 20 , whereas for 9 > |, there are only a finite 
number of persistent spins left. This important differ- 
ence is not adequately reflected in the pair correlation 
function, which has the same scaling form for all values 
of 9. In the next section, we study another quantity to 
characterize the spatial distribution which undergoes a 
rather significant change in its scaling properties across 
the transition. This quantity is the empty interval dis- 
tribution, which is one of the standard tools in the study 
of one-dimensional reaction-diffusion processes. 



the study of one-dimensional problems, and has been suc- 
cessfully applied to study spatial distribution of persis- 
tent spins in A + A — > model. Under the IIA, the 
relation between n(k,t) and P 2 (r,t) is 



P 2 (r,t) = n(r,t) + P(t)~ 



dxn(x,t)P 2 (r-x,t) (19) 



It is convenient to express this relation in terms of 
the Laplace transforms C(p,t) = C(r,t)e~ pr dr and 
n(p,t) — n(s,t)e~ ps ds, where C(r,t) was defined in 
the previous section. Under these transformations, Eq.19 
maybe expressed in the form 



n{jp, t) 



P(t)C(p,t) 
l + C(p,t) 



(20) 



From the scaling form for P 2 (r, t) given by Eq.15 , we 
find that 



C(p,t) = P(t)VtI(q,t) 



(21) 



where q = p\/i, and I(q,t) = J t °°i f(x)e qx dx The 
lower limit is put as instead of zero to take care of 
possible small argument divergence in the scaling func- 
tion. 

Let us first consider the case where 9 < \: In this 
case the scaling function f(x) is intcgrable, so we put 
the lower limit in the previous equation as zero. Using 
Eq.17, we find that 



I(q,t) = a M g M - 1 7(l - 29, qa) + V 9a 



< ~ 2 (22) 



where j(a,x) = J* e~ t t a ~ 1 dt is the incomplete 
Gamma function. After substituting in Eq.20 and Eq. 
21, and taking the t — > oo limit (keeping q fixed), we find 



B. The Empty Interval Distribution 

An empty interval, in our convention, is the separation 
between two consecutive persistent sites. The empty in- 
terval distribution (EID) n(k,t) is defined as the number 
of such intervals of length k at time t. For convenience, 
we also divide this quantity with the system size N so 
that n(k, t) satisfies the following normalization condi- 
tions. 



i; 



n(k,t)dk = P(t) 



kn(k,t)dk = l (18) 



Computing the EID directly, even under the mean field 
approximation, is non-trivial. Instead we shall compute 
it from the pair correlation function using the Indepen- 
dent Interval Approximation (IIA), where the lengths of 
successive empty intervals are considered as independent 
random variables. The IIA has been a valuable tool in 



n(p, t) = t i 



P(t)t* - 



a 2 V 8 7(l - 29,qa)+e-i a 



It follows that 



n(k,t) 



J 2 



so that 



Hp,t)=t-v[tvP{t)-h 1 (q)] 



(23) 



(24) 



(25) 



which has the same form as Eq.23, and h\(q) = 
J °° h(x)[l — e~ qx ]dx. After integrating by parts, we find 



hi(q) = -G(oo) + q lim xG(x) + G(q) 



(26) 



where G{x) = h{y)dy and G(q) = J °° G{x)e' qx dx. 
We assume that G(x) is integrable, so that 
lim x ^o zG^x) = and G(oo) = 0. Finally we have 
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G(q) 



-,26>„20 



q 2S j{l-26,qa)+. 



(27) 



Now we may try to deduce the behaviour of the func- 
tion G(x) at large and small arguments from its Laplace 
transform, given by the previous equation. To find the 
behaviour of G(x) near x — 0, we use the standard for- 
mula [18] 



lim t ^ t p g{t) = lim s 



9{s) 



;p>-l (28) 




where g(s) is the L.T of g(t). Now, for large q, 
7 (1 -26», qa) ~ r(l — 26>), so that G(q) = T7T=2gT as 1 ~* 
oo. It follows that lim x ^ x 1 ~ 2e G(x) 
from which we find 



(20-l)!r(l-20) 



G(x) 



(26 - l)W(l - 29) 



x->0 



(29) 



and, after using the relation h(x) = — |p, 



(30) 



The behaviour of G(x) at large x, one has to look for 
the singularities of G(q) in Eq. 27. If G(g) has a singular- 
ity of the form G(q) <~ (q — q*)~", then, upon inversion of 
the L.T, it follows that G(x) - x^e 9 ** as x -> oo [18]. 
In order to find the singularity, we plotted the denomi- 
nator of Eq. 27 against its argument (Fig. 2). We find 
that the function crosses zero at one point in the nega- 
tive q axis. By careful numerical analysis using bisection 
method, we have determined this crossing point to be at 
q*a = —A, where the numerical constant A ~ 0.32 for 
# = § and A ~ 0.85 for 6 = \. This implies that the 
leading term in the decay of G(x) at large x is exponen- 
tial, ie., G(x) ~ exp(— —x) as x — > oo, with a possible 
power-law prefactor. Consequently, the scaling function 
h(x) also has similar exponential decay at large x. 



FIG. 2. The figure shows the inverse of G(q) plotted against 



qa. 



We also determine the characteristic length scales of 
the distribution using the scaling form for n(k,t). These 
may be defined as the ratios of moments of the distribu- 
tion. 



L m (t) = I m (t)/I m -i(t) ;m = 0, 1,2, . 



(32) 



where I m (t) = ^2k n (k't)k m are the moments of 
n(k,t). Clearly, Li(t) = P(t)- 1 ~ t B by definition, 
whereas all higher order length scales 



Lj-(t)~t3 ;j>l 



(33) 



which follows from the dynamic scaling form given by 
Eq.24, 30 and 31 for n(k,t). 

We now continue our study of empty interval distri- 
bution for the case where 6 > \. Although the MMFA 
allows us to study arbitrarily large values of 9, we re- 
strict ourselves to the regime 6 < 1, since our basic aim 
is to study the persistence in g-state Potts model where 
9{q) < 1- Furthermore, for q = oo where 9 = 1, n(k,t) 
can be found exactly [7] and is known to be a pure expo- 
nential. 

For \ < 9 < 1, the scaling function f(x) has a 
non-integrable x~ 2e singularity near x = (We do not 
study explicitly the logarithmic singularity occurring for 
9 = 5). Wc integrate by parts and find 



,26» 



26 



-[<" 



,20-1 



7(2- 26, qa)] + -e~ qa ;6 > 



(34) 



h{x) 



~ e 



;x > 1 



(31) 



Let us now define A = -p^y = qt 2 p(t) 1 . For t — > 00 
and finite A, we have 



G 



I(q,t) = 



1 



29 -I ' \thp(t) 
After substitution in Eq.21 we find 

f 



C(p,t)=P(t) 



-,28 



29 - 1 



t e + 



XP(t) 



(35) 



(36) 



t e P(t)a 2e 

We define the constant (3 = — gF^i — i m terms of which 
C(p,t) =13+ j. Now we substitute in Eq.20 and find 



h{p,t) = P(t)- 



x + p- 1 



•1 + 

Upon inversion of the L.T, we find that 
n(k,t) = P(t) 2 (j)[kP(t)} 

where 



4>{x) 



1 + 



8(x) + 



1 



(3{l + (3) 



-x(l+/3)- 



(37) 



(38) 



(39) 



The scaling function has a rather unnatural <5-function 
singularity at the origin. However, a more careful anal- 
ysis show that for any finite (but still large) time t, the 
divergence at origin is only power-law, but with a differ- 
ent exponent than the previous case (9 < i). We start 
with the expression given by Eq.21 and Eq.34 for C(p, t). 
After keeping the leading finite t correction, we find that 

c(p,t) = p + \-pr e ^ e - 1 h 2e - 1 (40) 

A 

For purposes that will be clear later, let us define 
m(k 1 t) = kn(k, t) so that m(k 1 1) — 1. We also define 
the associated Laplace transform rh(p,t). The Laplace 
transforms are related through 



Mp, t) = - 



dn{p, t) 
dp 



(41) 



Using the expression Eq.20 for h(p,t), we find that 
P(t)C'(p,t) 



m{p, t) 



where C'{p,t) - 



(42) 

[i + c{p,t)f 

'j^ 1 , and is given by the expression 



C'(p,t) 



P(t) 



^ + (29-l)l3r e ^-^\ 2 ^ 
A 



(43) 



After substitution in Eq.41 and taking the limit t 
oo, we find 



fhijp, t) 



l + (29-\)(3t-^ 2e - 1 )\ 20 
[l + A(l+/3)] 2 



(44) 



which gives the scaling forms 

m(k, t) = P(t)ip(t, kP(t)) ; n(k, t) = P(t) 2 $(i, kP(t)) 

(45) 

where 



x$(t, x) = i/)(t,x) 



(46) 



by definition. The Laplace transform of the scaling 
function ip(t, x) is 



#,A) = 



l + (29-l)f3t- 9( - 2e ^\ 2e 
[1 + A(1 + /?)] 2 



(47) 



We notice that if the finite t correction term is not in- 
cluded, limx^oaX 2 ^^, A) is finite, and in that case, the 
small argument divergence of tp(t, x) will be sharper than 
any power-law. This is what is reflected in the appear- 
ance of the 5-function in Eq.41. However, when this term 
is included, the multiplying factor has to be A 2-29 in or- 
der to make the resulting expression finite as A — > oo. 
This implies that the small x divergence for ij){t,x) has 
the power-law form ip(t, x) <~ i-6(2B-i) x i-2e f or sma u x 
From Eq. 48, we find a similar power-law divergence in 
<&(t, x) also. 



Ht,x) 



t 



-e(20-i) x -20 



\x < 1 



(48) 



In the large x limit, $(i, x) becomes time independent, 
and decays exponentially with x as in Eq. 39, ie., 



Ht,x) 



(1 + /3) 2 



e !+" 



;x > 1 



(49) 



The characteristic length scales are easy to compute. 
From the scaling form, it follows that all the characteris- 
tic lengths have identical asymptotic scaling behaviour. 



Lj(t)~t ;i = 1,2, 



(50) 



The difference in the asymptotic scaling behaviour of 
the characteristic length scale as 9 crosses \ may be seen 
as a competition between two length scales, the diffusive 
length scale Co(t) ~ \[Dt which gives the mean sepa- 
ration between two random walkers, and the persistence 
scale C p {t) = Pit)^ 1 ~ t e which is the mean separation 
between two persistent spins. The characteristic length 
scale is dominated by the larger of the two, ie., we may 
write 



=max(i 6) 



(51) 



where C(t) is defined through the dynamical scaling 
form for n(k, t). 



n(k,t)=£(t)- 2 g(t,-^) 



(52) 
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The scaling function g(t,x) = h(x) when 9 < 5 and 
g{t,x) — &(t,x) when 9 > \. In general, the small argu- 
ment behaviour of g(t, x) has the power-law form 

g(t,x) ~ r^x- T ;x^0 (53) 

where the exponents ip and r are given by 

^ = 0(20-l)JT(0-^) ; r =max[26>, 2(1 - 0)], (54) 

and i?(x) is the Heaviside step function. For large x, 
the scaling function is time-independent and decays ex- 
ponentially with x. We also note from the scaling form 
that over small distances k <i*, 

n{k,t) <~ P(t)k- T ;fc<Ct (55) 

where r is given by Eq.54. 

IV. NUMERICAL RESULTS 

We studied the quantities P2(r,t) and n(k,t) numeri- 
cally by simulating the kinetics of g-state Potts model 
with random initial conditions. The time evolution 
of spin configurations via Glauber dynamics is imple- 
mented using the mapping of this dynamics to a reaction- 
diffusion problem, as mentioned in the introduction. In 
this procedure, a set of diffusing particles A are ini- 
tially distributed at random on the lattice with a cer- 
tain average initial density no (which we fix as 5). When 
two diffusing particles meet, they annihilate each other 
or coagulate with probability ^-j- and respectively. 
We count one MC step in the simulation after the po- 
sition of every particle in the lattice has been updated 
once. Persistent spins(sites) at any time t are those sites 
which have not been touched by a random walker till 
that time. All the simulations were done on a lattice 
with 2 17 sites, and the results were averaged over 100 
different starting configurations. In order to check the 
different dynamic scaling behaviour for 9 < 5 and 9 > 5 , 
we did our simulations for three different values of q- 
2,5 and 10. For later reference, we note that from Eq.l, 
the corresponding values of the persistence exponent are 
6>(2) = 3/8 = 0.375,0(5) w 0.6928 and 9(10) w 0.8310. 
In Figs. 3-5, and later in Figs. 7-10, we have employed 
logarithmic binning of the data in intervals of size 1.5™ 
(n = 1,2,....). since the statistical noise was consider- 
able. However, for all exponent measurements, we have 
used only the bare (not binned) data. 

In Fig. 3-5, we have plotted the scaling function f(x) 
for the pair correlation function P 2 (r, t) against the scal- 
ing variable rj\ft for three q- values, q — 2,5 and 10. We 
find excellent scaling collapse for all three values of q, 
which is in agreement with the dynamic scaling picture 
provided by the MMFA in Eq.15. In the figures, we find 



power-law decay of f(x) for small x, with a sharp cross- 
over to the flat long-distance behaviour, which is also in 
agreement with the assumption we made in Eq.17. We 
also note that the constant a introduced in Eq.17 is in 
fact very close to 1. 



. 1 II 1 II 


1 II II II II. 

t = 10 3 — : 




t = 8 x 10 3 




t = 32xl0 3 ■ 




\ slope = -0.75 ■ 


, , ,1 , , ,1 


, , ,1 , , ;>. , , ,1 , , , 



10- 3 10- 2 10- 1 1 10 10 2 10 3 

1 = rj\ft 



FIG. 3. The scaled pair correlation f(x) = P{t)~ 2 P 2 {r,t) 
is plotted against the dimensionless scaling variable x — r j\ft 
for q — 2 on a logarithmic scale. The straight line is a guide 
to eye, and has slope 261(2) = 0.75 , which is the MMFA pre- 
diction. The time t is measured in MC steps and distance r 
is measured in units of lattice spacing. 





IJ 1 1 II 1 
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1 l| 1 1 l| 1 II; 

t = 10 3 — : 




' \ 




t = 8xl0 3 , 








t = 32 x 10 3 








slope = -1.38 




,1 , , ,1 , 


\ \^=p^ — 

, ,1 ,\ , ,1 


, ,1 , , ,1 , , , 



10- 3 10- 2 10- 1 1 10 10 2 10 3 10 4 

x = r/yfl 



FIG. 4. Same as Fig. 3, for q — 5. The slope of the straight 
line is 20(5) ~ 1.38, which is the MMFA prediction. 
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FIG. 5. Same as Fig.3, for q = 10. The slope of the straight 
line is 20(10) ~ 1.67, which is the MMFA prediction. 

In Fig. 6, we plot the characteristic length scale L 2 (t) 
against time t for q = 2, 5 and 10 and measure the dy- 
namical exponent <p. The observed slopes of the lines 
are systematically higher than the theoretical prediction 
in Eq.53 by around 0.05, while the statistical error in 
all the three cases was only <~ 10~ 4 or smaller. The ob- 
served deviation could be possibly due to the fact that the 
asymptotic behaviour is not fully reflected over the time 
scales which we used. The number of persistent spins left 
in the system falls rapidly with time for high values of 
q, and so we were forced to restrict ourselves to times 
t < 32000. In fact, even for q = 2 case, previous simula- 
tions over longer time scales have shown the presence of 
an additive power-law correction to the asymptotic scal- 
ing behaviour [6] . 




FIG. 6. The figure shows the characteristic length scale 
Z/2 (i) (measured in units of lattice spacing, definition in text) 
plotted against time t (measured as number of MC steps) on a 
logarithmic scale for three Potts values q = 2, 5 and 10. The 
measured slopes of the lines are respectively 0.5507, 0.7391 
and 0.8672. The corresponding theoretical predictions are, to 
the same accuracy, 0.5000, 0.6928 and 0.8310. 



FIG. 7. The scaled EID h(x) = tn(k, t) is plotted against 
the dimensionless scaled separation x = k/\ft for q = 2 
Potts model. The excellent scaling collapse validates the 
scaling form given in Eq. 24. The straight line has slope 
r = 2(1-61(2)) = 5/4, which is the MMFA-IIA prediction.The 
time t is measured in MC steps and distance k is measured in 
units of lattice spacing. 



In Fig. 7-9, we check the dynamic scaling form Eq.54 
for n(k,t) against the scaling variable x = h/fr? for three 
values of 9-2,5 and 10. We find that for q = 2, excellent 
scaling collapse is obtained with (j) — \ (Fig. 7). For 
small x, we find power-law decay of the scaling function, 
which crosses over to fast exponential decay at large x. 
For higher values of q (where 6{q) > |), we find that 
with <p = 9, we find very good scaling collapse for x S> 1. 
But for x <C 1, we find systematic deviation from scal- 
ing collapse, which was also observed earlier by Bray and 
O'Donoghuc [7]. This observation supports the theoret- 
ical prediction based on MMFA, and shows that in this 
regime, the scaling function has explicit time dependence. 
To show this more clearly, and to verify the predicted 
time-dependence, we plotted the quantity n(k,t)/P(t) 
against k for three widely spaced values of t for all three 
q values studied in Fig. 9. We see that in all three cases, 
a simple power- law decay with k is observed for k 
thus validating Eq.55. The measurement of the exponent 
r gives values in reasonable agreement with theoretical 
prediction, although for q = 10, the statistical error is 
significant. 
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FIG. 8. The scaled EID x) = tn(k, t) is plotted against 
the scaled separation x = kP(t) on a logarithmic scale for 
q — 5 Potts model. We see that the scaling function is explic- 
itly time-dependent for small k, but is time- independent for 
large k. The straight line in the figure gives the theoretical 
prediction r = 29(5) ~ 1.38 for power-law decay at small x 
(see discussion in text). The time t is measured in MC steps 
and distance k is measured in units of lattice spacing. 



FIG. 10. In the figure, n(k,t) / P(t) is plotted against k for 
two widely separated values of t for each value of q = 2, 5 and 
10 (top to bottom). In all the cases, the function is indepen- 
dent of t for k <JC i*, and shows the power law decay ~ kT T . 
We measure r ~ 1.32 ± 0.03, 1.41 ± 0.04 and 1.61 ± 0.22 for 
q — 2, q — 5 and q = 10 respectively. The corresponding 
MMFA-IIA predictions are, to the same accuracy, 1.25, 1.38 
and 1.66. 



: 1 ' V 1 




1 1| i 1 1| i ii. 
t = io J — - 






f = 8x 10 3 ■ 




A 1 


= 32 x io 3 — ; 
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FIG. 9. Same as Fig. 6, for q = 10. The straight line in 
the figure gives the theoretical prediction r = 26(10) ~ 1.67 
for power-law decay at small x (see discussion in text). The 
time t is measured in MC steps and distance k is measured in 
units of lattice spacing. 



V. CONCLUSIONS 

In this paper, we have studied the spatial aspects of 
persistence in one dimensional (/-state Potts model using 
a mean field approximation. We have computed the pair 
correlation function for persistent spins under this ap- 
proximation, and used it to compute the empty interval 
distribution under the independent interval approxima- 
tion. We find dynamical scaling behaviour in both these 
quantities. The time dependence of the characteristic 
length scale and the behaviour of the scaling function was 
found in both cases. We showed analytically within the 
mean field approximation the transition from fractal to 
homogeneous distribution of persistent spins as the per- 
sistence exponent crosses |. We support our results by 
numerical simulations in the kinetic g-state Potts model. 
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